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We search for pair-produced Dirac magnetic monopoles in 35.7 pb _1 of proton-antiproton col- 
lisions at ^/s = 1.96 TeV with the Collider Detector at Fermilab (CDF). We find no monopole 
candidates corresponding to a 95% confidence-level cross-section limit a < 0.2 pb for a monopole 
with mass between 200 and 700 GeV/c 2 . Assuming a Drell-Yan pair production mechanism, we set 
a mass limit m > 360 GeV/c 2 . 

PACS numbers: 



The existence of magnetic monopoles would add sym- 
metry to Maxwell's equations without breaking any 
known physical law. More dramatically, it would make 
charge quantization a consequence of angular momentum 
quantization, as first shown by Dirac With such ap- 
peal, monopoles continue to excite interest and have been 
the subject of numerous experimental searches. 

Grand unified theories predict monopole masses of 
about 10 17 GeV/c 2 , so cosmic ray experiments have 
searched extensively for high-mass monopoles produced 



in the early universe. Accelerator searches for low-mass 
monopoles have looked for the effects of virtual monopole 
loops J3, EI Q 1 but the results have been questioned [j| . 
Detector materials exposed to radiation from pp colli- 
sions at the Tevatron have been examined for trapped 
monopoles, but the limit obtained depends on the model 
for the trapping of monopoles in matter |(| . Despite these 
efforts, magnetic monopoles have not been discovered Q. 

Magnetic monopoles have magnetic charge g satisfying 
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the Dirac quantization condition: 

— = - - = — 68 5 • n 

he 2 e 2a 

where n is an integer and a is the fine structure constant. 
In this search, we consider an n = 1 monopole with mass 
less than 1 TeV/c 2 , spin A, and no hadronic interactions. 
Monopoles are accelerated by a magnetic field and are 
highly ionizing due to the large value of g/e. 

This search uses a 35.7 pb _1 sample of pp collisions at 
\/s = 1.96 TeV produced by the Fermilab Tevatron and 
collected by the CDF II detector during 2003 using a 
special trigger. The detector consists of a magnetic spec- 
trometer including silicon strip and drift-chamber track- 
ing detectors and a scintillator time-of-flight system, sur- 
rounded by electromagnetic and hadronic calorimeters 
and muon detectors CDF uses a superconducting 
solenoid to produce a 1.4 T magnetic field. The field is 
parallel to the beam direction, which is taken as the z 
direction, with c/> the azimuthal angle, and r the radial 
distance in the transverse plane. 

The important detector components for this search are 
the central outer tracker (COT) and the time-of-flight 
(TOF) detector E3 > both positioned inside the solenoid. 
The coverage of the cylindrical COT extends from a ra- 
dius of 40 cm to 137 cm and to psuedo-rapidity \rj\ ~ 1. 
The COT consists of eight superlayers, each containing 
12 layers of sense wires. The COT makes timing mea- 
surements for track reconstruction as well as integrated 
charge measurements for determining a particle's ioniza- 
tion energy loss dE/dx. The COT is surrounded by 216 
TOF scintillator bars, which run parallel to the beam 
line and form a cylinder of radius 140 cm. Each TOF 
bar is instrumented with a photomultiplier tube (PMT) 
on each end. The TOF measures both the time and 
height of PMT pulses; the pulse height is typically used 
to correct for discriminator-threshold time-slewing. Due 
to their large ionization and massive production of delta 
rays, monopoles in scintillator with velocity f3 > 0.2 are 
expected to produce more than 500 times the light from 
a mimmum-iomzmg particle (MIP) 0E3- 

We have built and commissioned a highly ionizing par- 
ticle trigger that requires large light pulses at both ends 
of a TOF scintillator bar. The trigger was designed to de- 
tect monopoles efficiently while consuming less than 1 Hz 
of the CDF data acquisition bandwidth. The electron- 
ics response of the TOF has been calibrated EE3 to 
account for non-linearities and channcl-to-channcl differ- 
ences. The trigger thresholds of about 30 MIPs are well 
below the expected response to a monopole and have a 
negligible effect on the trigger efficiency. 

In the CDF detector, a monopole is accelerated along 
the uniform solenoidal magnetic field in a parabola 
slightly distorted by relativistic effects. Because no other 
particle mimics this behavior, the TOF acceptance must 
be estimated from Monte Carlo simulation. We have ex- 



tended the GEANT simulation [II El El [13 to handle 
magnetic monopoles, including the acceleration from the 
magnetic field, energy loss and multiple scattering |l7j . 
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FIG. 1: The acceptance of the TOF for monopole pairs, as 
a function of monopole mass. The band indicates the total 
systematic uncertainty. 



Because the monopole-photon coupling is large and 
non-perturbative, there is no universally accepted field- 
theoretic calculation of magnetic-monopole production. 
However, monopole interactions with matter, such as 
scattering, require only a replacement of the electric 
charge with the monopolc's effective charge g(3. This has 
led the authors of Ref. j(| to adopt a heuristic production 
model by making the same replacement for Drell-Yan 
monopole pair production, which we take as our primary 
benchmark. 

Either a monopole or anti-monopole must reach the 
TOF detector in order to cause a trigger. To calcu- 
late the TOF acceptance for the heuristic pair produc- 
tion mechanism, we produce lepton Drell-Yan events 
with Pythia E3 with the lepton mass replaced by the 
monopole mass, and weight events according to the ad- 
ditional velocity dependence. The TOF acceptance for 
monopole pairs simulated with GEANT is shown in Fig- 
ure n] Light monopoles, accelerated strongly by the 
magnetic field, tend to be swept out of the detector be- 
fore reaching the TOF. Heavy monopoles, produced near 
threshold, suffer the same fate. 

Because we are unable to test experimentally the 
model for material interactions, we assign a systematic 
uncertainty of one half the total calculated effect, mea- 
sured by comparing the TOF acceptance for the full sim- 
ulation with a ficticious detector consisting of the TOF 
only. The material in the detector lowers the acceptance 
due to energy loss and multiple scattering, a 3% system- 
atic error for intermediate-mass monopoles. This method 
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likely overestimates the uncertainty; varying the energy- 
loss model between a naive model where e — > and the 
full treatment of Ref. has a negligible effect. 

The TOF acceptance depends on the monopole pro- 
duction kinematics. To quantify this dependence, we 
consider separately the Drell-Yan mechanism without the 
additional velocity dependence and with monopole pro- 
duction uniform in the cosine of the polar angle in the 
center of mass frame. The total variation in the accep- 
tance is 10%. We therefore present results for our bench- 
mark mechanism only, with the understanding that mass 
limits for other production mechanisms can be inferred 
from the cross-section limit with reasonable accuracy. 

During each event, the TOF electronics makes a sin- 
gle measurement for each PMT. Light from other par- 
ticles, called spoilers, can reach a PMT before the light 
from monopoles, starting the charge integration. If the 
monopole light does not reach the PMT within the 20 ns 
charge integration window, the monopole's light will not 
be integrated and trigger will not fire. Our studies show 
that pure Monte Carlo underestimates the effect of spoil- 
ers seen in data. We therefore estimate the spoiler frac- 
tion by embedding Monte Carlo produced monopoles in 
real Z — > e + e~ data. Because these are high- mass cen- 
tral events produced by a Drell-Yan mechanism, we ex- 
pect the distribution of other particles in the event to 
be similar to that of a monopole-pair production event. 
We exclude the bars with signals from the electrons and 
count the number of spoiler events, which have real pulses 
arriving more than 20 ns before the simulated pulse from 
a magnetic monopole. 

The systematic uncertainty is dominated by the un- 
certainty in the time needed to integrate enough of the 
monopole's charge to cause a trigger. To quantify this 
effect, we note that rise times for TOF pulses are typi- 
cally less than 1 ns and redo the calculation with a 15 ns 
integration window. We take one-half the difference as 
a systematic uncertainty. Other effects, such as the de- 
pendence on luminosity, are much smaller for our sam- 
ple. For a 400 GeV/c 2 monopole, the spoiler fraction is 
2% ± 1% with a 3% systematic uncertainty. 

Massive monopoles can have low velocities causing 
them to arrive at the TOF too late to cause a trigger. The 
timing acceptance is calculated with a Monte Carlo sim- 
ulation by requiring pulses to arrive within the 54 ns tim- 
ing window. Only heavy monopoles move slowly enough 
to be affected: a 900 GeV/c 2 monopole is out of time 
in 10% of events. This is a negligible effect on lighter 
monopoles. 

Monopoles curve in the rz plane, in sharp contrast to 
electrically charged particles, which curve in the r<p plane. 
A specialized reconstruction program isolates monopole 
candidates using data from the COT. Candidates consist 
of coincident track segments composed entirely of hits 
with large ionization, consistent with a straight line in 
the plane perpendicular to the magnetic field. 



The COT electronics encodes the integrated charge as 
the width of a hit, which is the ionization measurement 
used for monopole candidate selection. A typical MIP 
produces hit widths of about 20 ns. An extrapolation of 
the non-linear COT response for ordinary particles pre- 
dicts that monopoles would produce hit widths of about 
230 ns (1000 MIPs), still within the dynamic range of the 
COT. We do not use this extrapolation. Instead we cut 
in the tail of the width distribution from ordinary tracks, 
found to be at 140 ns (50 MIPs) in minimum-bias data 
collected with an open trigger highly efficient for inelastic 
pp collisions. Hits with charge below this amount are not 
considered by the monopole reconstruction. As magnetic 
monopoles have much greater ionization than the tracks 
used to determine this cut, it has a negligible effect on 
the efficiency. 

The default COT tracking algorithm first reconstructs 
track segments in each of eight superlayers. It checks for 
hits loosely consistent with a straight-line, using a tol- 
erance of 20 ns. The identified hits in each segment are 
then fit to a circular trajectory. In the monopole algo- 
rithm, the segments are required to be composed entirely 
of high-ionization hits. Also, because a monopole can be 
as slow as (3 ~ 0.1 with changing transverse velocity, the 
usual timing assumption (iflight = r/c) cannot be used. 
Instead, the time of flight to each superlayer is varied 
between r/c and lOr/c in 5 ns increments. 

A monopole candidate consists of several 0-coincident, 
low-curvature segments. From Monte Carlo simulation, 
we choose a loose cut on the segment curvature p < 
0.001 cm -1 , which for an electron would correspond to 
Pt > 4 GeV/c. Likewise, the <fi tolerance is a loose 
0.2 radians. The remaining cuts are on the minimum 
number of hits needed in a segment and on the total num- 
ber of (/(-coincident segments required for a monopole can- 
didate. By ignoring the width cut, the segment-finding 
algorithm efficiency is measured in an independent data 
sample using high-j»T tracks. In this manner, we choose a 
highly efficient cut requiring seven coincident superlayers 
with at least eight hits in each segment. This has a 94% 
efficiency with a 1% statistical uncertainty. For these 
cuts, the efficiency for finding high-mass monopole pairs 
calculated with the Monte Carlo simulation is nearly 
100%. The efficiency for high-px electrons in simulation, 
after removing the width cut, is also nearly 100%. There 
are real detector effects contributing a small inefficiency. 

As an ionizing particle passes through matter, the most 
energetic electrons form delta rays. For highly relativis- 
tic low-mass monopoles, the large number of delta rays 
confuses the segment finding algorithm, lowering the ef- 
ficiency. We check that GEANT is properly producing 
delta-rays by comparing the efficiency of monopoles to 
kinematically equivalent heavy-ions simulated in the ab- 
sence of a magnetic field. We scale the efficiency deter- 
mined from Monte Carlo simulation to make the high- 
mass monopole efficiency agree with the high-px track 
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FIG. 2: Number of monopole candidates as a function of 
COT width cut in offline reconstruction, in 100k minimum- 
bias events and the entire 130k event trigger sample. A 100 
ns width cut corresponds to 15 MIPs, we expect monopoles 
to ionize more than 1000 MIPs or 232 ns. 



FIG. 3: The 95% CL cross-section upper limit versus 
magnetic monopole mass. The theory curve for Drell-Yan 
monopole pair production intersects at the mass limit m > 
360 GeV/c 2 . 



Effect 


Efficiency 


TOF geometric (MC) 
TOF response 
TOF spoilers 
TOF timing (MC) 


70% ± 3% ± 3% 
100% 

98% ± 1% ± 3% 
99% ± 1% ± 1% 


COT width cut 
COT segment finding 


100% 

94% ± 1% ± 3% 



TABLE I: Efficiency of the monopole search with statis- 
tical and systematic uncertainties for a monopole mass of 
400 GeV/c 2 . The full mass dependence is accounted for in 
the limit. 



efficiency. As the small inefficiency from real detector ef- 
fects cannot be measured directly on monopoles, we take 
one half of the total inefficiency as a systematic uncer- 
tainty: 3% for 400 GeV/c 2 monopoles. 

To estimate how effectively the monopole reconstruc- 
tion rejects background, we use minimum-bias data. In 
8 x 10 5 events, the event most like a monopole has two 
coincident super-layers with seven hits per segment. Our 
monopole requirements are much more stringent. We 
require a seven-fold coincidence of eight hits or more, re- 
sulting in extremely small background. In the trigger 
sample the background is similarly small; the event most 
like a monopole has two coincident super-layers with six 
hits per segment. In Fig. we count the number of 
monopole candidates passing looser cuts on the hit width. 

None of the 130, 000 events from the monopole trigger 
sample passes the candidate requirements, and we report 
a limit |l3( . Monopole production limits are typically re- 
ported by the cross-section upper limit as a function of 



monopole mass to minimize the dependence on a partic- 
ular production model. The expected number of events 
N from a process with cross section a and detector effi- 
ciency with acceptance e after integrated luminosity L is 
given by N = Lea. We calculate the cross-section limit 
for zero observed events, based on the efficiency summa- 
rized in Table [Hand a 6% uncertainty in the luminosity 
measurement |l9j . We find the cross section for which 
pseudo-experiments with efficiency and luminosity cho- 
sen randomly according to their uncertainties yield one 
or more measured events 95% of the time. 

Our cross-section exclusion limit is shown in Figure 
Our limit excludes monopole pair production for cross 
sections greater than 0.2 pb at the 95% confidence level 
for monopole masses between 200 and 700 GeV/c 2 . For 
the Drell-Yan mechanism, this implies a mass limit of 
to > 360 GeV/c 2 at the 95% confidence level. This 
is currently the best limit from a direct search. Addi- 
tional Run II data will improve the sensitivity: another 
300 pb -1 extends the mass reach by 100 GeV/c 2 . 
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